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Abstract 

We discuss the dynamics of two harmonic oscillators of which one has a negative kinetic 
term. This model mimics the Hamiltonian in quantum geometrodynamics, which possesses 
an indefinite kinetic term. We solve for the time evolution in both the uncoupled and coupled 
case. We use this setting as a toy model for studying some possible aspects of the final 
stage of black-hole evaporation. We assume that one oscillator mimics the black hole, while 
the other mimics Hawking radiation. In the uncoupled case, the negative term leads to a 
squeezing of the quantum state, while in the coupled case, which includes back reaction, we 
get a strong entangled state between the mimicked black hole and the radiation. We discuss 
the meaning of this state. We end by analyzing the limits of this model and its relation to 
more fundamental approaches. 
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1 Introduction 



In spite of many attempts, a final theory of quantum gravity remains elusive [T]. Such a theory is. 
however, needed for at least two reasons. First, there are conceptual and formal arguments which 
point to an encompassing fundamental framework. And second, such a theory is required in order 
to tackle concrete physical problems. Among the most important ones are problems in cosmology 
and black-hole physics. In our paper we address the latter field. There, two main issues are usually 
studied in this context: the microscopic interpretation of black- hole entropy and the description 
of the final evaporation phase. The first problem is already relevant for large black holes, that 
is, for black holes with a mass much bigger than the Planck mass. Therefore, some progress has 
been achieved in various approaches such as string theory, loop quantum gravity, and quantum 
geomctrodynamics pQ[21[3]. On the other hand, as the black hole approaches the Planck mass, a 
full understanding of the theory is necessary for the final evaporation phase, and it is therefore 
not surprising that no definitive conclusion has been reached up to now. 

In our paper we attempt to shed some light on the question of how the final phase in the 
black-hole evolution may be described in terms of quantum states used in canonical quantum 
gravity. In this approach, the central kinematical entities are states that depend - in the gravi- 
tational sector - on the three- metric (in quantum geometrodynamics), a non-abelian connection 
(in quantum connection dynamics), or a SU(2)-holonomy (in loop quantum gravity) The 
central 'dynamical' equations are the Hamiltonian and the diffeomorphism constraints (plus, in 
the latter two versions, the Gauss constraint). The final stage of black- hole evaporation has not 
yet been described in this language. There exist several attempts using the semiclassical Einstein 
equations, G^i, — 87rG(T^^), in order to implement the back reaction of the Hawking radiation on 
the evaporating black hole, with no final result. Since the semiclassical Einstein equations cannot 
be fundamentally correct [l] , a definite answer can only be obtained from the exact theory. It is the 
purpose of this paper to provide a possible insight into the quantum states of the last evaporation 
stage. We are certainly not able to present a realistic description of the final stage. But we can 
address and tentatively answer the following question: Suppose we have solved for the quantum 
state - how can we recognize black hole evaporation from it? We achieve this goal through some 
admittedly oversimplified models which, however, as we hope, would capture some real features 
of the physical process under investigation. The model is mainly chosen in order to study the 
influence of an indefinite kinetic term, as it occurs in the Wheeler-DcWitt equation. Thus, in this 
respect our model may exhibit realistic features of the black-hole evaporation process, while in 
other respects it may be totally unrealistic and must be replaced with features from more realistic 
(at the moment not exactly soluble) models. We herein envisage the black hole as being embedded 
as a quantum object into a semiclassical universe. We thus have to address the time-dependent 
Schrodinger equation with the full quantum gravitational black hole Hamiltonian mimicking the 
black hole. 

Our paper is organized as follows. Section 2 presents a brief review of the semiclassical expan- 
sion. We present our arguments why we use a (functional) Schrodinger equation with a quantum 
gravitational Hamiltonian aiming to describe a quantum black hole. Section 3 is the main part 
of our paper. We present our model of coupled oscillators mimicking the quantum black hole 
and Hawking radiation. We first consider the case without direct back reaction (only indirectly 
through the decrease of the black- hole mass). The negative kinetic term (characteristic of a black- 
hole Hamiltonian) leads to a squeezed state when the mass which we associate with a black hole 
approaches the Planck mass. We then discuss the case of a direct coupling. This leads to an 
entangled state between the mimicked black hole and the Hawking radiation. We calculate and 
discuss the corresponding reduced density matrices. Section 4 contains our conclusions. 

2 Semiclassical limit and beyond 

The central equations of canonical quantum gravity in the geometrodynamical, connection, or loop 
approach are of the constraint form iJ^P = 0. We shall restrict in our paper attention to quantum 
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geometrodynamics; from loop quantum gravity one would expect in addition some features arising 
from the discreteness of space, which are not captured here, but which could play a crucial role in 
a more realistic scenario [5]. The first point to notice is that the constraint equations in the form 
77 \1/ = are not yet suitable to describe black-hole evaporation. This is because these equations, 
as they stand, describe a closed system where everything is described quantum gravitationally. In 
order words, they would describe the case of a quantum black hole within a quantum universe. 
Although this may be the appropriate picture at the most fundamental level, the situation that 
one wants to address is a quantum black hole embedded in a semiclassical universe for which an 
appropriate time parameter is present. 

Such a time variable can arise in various ways. One can, for example, discuss the gravitational 
collapse of a dust cloud, in which a dust proper time emerges in a natural way, cf. [1] and the 
references therein. Solutions of quantum geometrodynamics can be exploited to derive Hawking 
radiation and greybody corrections within this model [5j. Another way of recovering time is a 
Born-Oppenheimer type of expansion scheme at the full (formal) level of quantum geometrody- 
namics [T]. This is the framework which we shall use here. 

Performing, for example, an expansion with respect to the Planck mass, one can derive from 
the full constraint equation (the Wheeler-DeWitt equation) a functional Schrodingcr equation for 
'matter' fields in an external spacetime. The time in this Schrodingcr equation is of a semiclassical 
nature and defined by configurations of the semiclassical gravitational field. More concretely, this 
equation reads 

I d'x {a^(x)7i:T(x) + A^'^(x)K"(x)} . (1) 

Usually, iJ™ is the Hamiltonian for the 'matter' (i.e. non-gravitational) fields in the Schrodinger 
picture, parametrically depending on (generally non-static) metric coefficients of the curved space- 
time background. The 'bra' and the 'kef in the quantum state refer to the Hilbert space of the 
'matter' fields. {N and N'^ are lapse and shift function, respectively, and 'H™{x.) and 7i™(x) denote 
the quantum 'matter' part of the Hamiltonian constraint and diffeomorphism constraint operator, 
respectively.) It is important to emphasize that only the semiclassical degrees of freedom of the 
gravitational field enter the definition of t. 

In realistic situations, the full set of these semiclassical variables may, however, include only 
part of the gravitational field and may also include part of the matter. On the other hand, the 
gravitons, which are small excitations of the metric, behave fully quantum and must be included 
into i?™. But this is also what happens in our case here: The degrees of freedom corresponding 
to the quantum black hole must be part of while the time parameter t is defined by the 
macroscopic part of gravity and matter, that is, by the semiclassical Universe into which the 
quantum black hole is embedded. 

To describe black-hole evaporation, therefore, one should stick to this semiclassical description 
of the embedding universe - providing the time parameter t -, while employing a full quantum 
description for the black hole. This corresponds to the realistic situation of an observer residing 
outside the quantum black hole and having the usual semiclassical time at his disposal. In such 
a situation one has therefore to apply Equation ([T]) with t referring to the semiclassical time of 
the outside universe, and if™ being the full Hamilton operator of the quantum black hole and the 
fields interacting with it. It is for this reason that we consider an indefinite kinetic term in H™; 
it is inherited from the Wheclcr-DeWitt equation and applicable here because we apply quantum 
gravity to the black hole. One thus has to employ a mixture of Schrodinger and Wheeler-DeWitt 
equation0 

The kinetic term of the gravitational part of the quantum Hamilton operator is suppressed 
by the Planck mass, mp. As long as the black- hole mass is large, this kinetic term should thus 

^The quantum formation of a black hole from spherical domain- wall collapse was studied recently in a related 
framework in [6] . 
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be irrelevant. One has in this hniit the Hawking radiation as the only contribution to ([T]). The 
corresponding quantum state was explicitly calculated and discussed in [7j. After the black hole 
approaches the Planck mass in the final evaporation phase, the full quantum Hamiltonian of the 
black hole becomes relevant. In the full field theoretic framework of [7], this was not considered. 
We shall instead investigate this question here in a quantum mechanical model that can capture 
some of the relevant features. In particular, as emphasized above, we have to deal with the 
important property of the gravitational Hamiltonian possessing an indefinite kinetic term. This 
is important for the understanding of time in quantum gravity [H [8] , but will also be crucial for 
the qualitative features of black-hole evaporation. 

3 A simple model of black-hole evaporation 

The full equation ([1]) is a complicated functional differential equation for a wave functional de- 
pending on the three-metric and matter fields (including Hawking radiation). In order to make 
our analysis tractable, we shall instead consider the following model, which is purely quantum 
mechanical: 



How do we interpret this ansatz? Let us first recall that the time t in this equation is the 
semiclassical time of ([T]) coming from the scmiclassical degrees of freedom (such as the scale factor 
and macroscopic matter) of the Universe. Equation ((T|) contains second functional derivatives 
with respect to the three-metric describing the black hole. Wc shall mimic this three-metric by 
the single variable x; we expect x to be the mass, M, of a Schwarzschild black hole (more precisely, 
its Schwarzschild radius 2GM / c^). The 'wrong' sign of its kinetic term refiects this correspondence 
to the three-metric. At a stage where the black hole is quite large, we expect that the kinetic term 
referring to x is negligible. However, for the final phase which we shall discuss in this paper, M 
will be small. Therefore, such a term will be of relevance, which is why it has been included into 



The variable y is supposed to be the analogue to Hawking radiation, with my corresponding to 
its energy, and z stands for the remaining degrees of freedom (with to^ as a formal mass parameter) . 
We shall consider in the following only the x- and y-variables, that is, wc shall restrict attention 
to the black hole and its interaction with Hawking radiation. Wc assume that at least for large 
black holes the evolution is stable in the sense that the potential docs not become too negative. 
For simplicity we have chosen harmonic oscillator potentials; this is realistic for the description 
of the Hawking radiation [7], but definitely oversimplified for the black hole. In ^ we have not 
yet taken into account the direct back reaction of Hawking radiation (y-part) onto the black hole 
(cc-part), that is, no x-y-coupling is included in This back reaction will be implemented in 
Section 3.3 below. 

We admit that the whole ansatz is not fully realistic, but we hope that our results will shed 
at least some light on the interpretation of black-hole evaporation in quantum gravity: the set- 
ting associated with (2) is not that of a standard quantum mechanical description of harmonic 
oscillators. On the one hand, there is an unusual kinetic term. As we shall see, the main conse- 
quences in our analogue to black-hole evaporation come from the negative kinetic term for x and 
arc independent of the details of the potential. On the other hand, the Planck mass appears as 
a parameter (suggested from the full Whceler-DeWitt equation, where it appears in front of the 
derivatives with respect to the three- metric), participating in (2) as a natural 'suppressor'; this is 
different from the dynamics of the Hawking radiation [7] . 

Let us be more precise, taking only the x- and y-degrees of freedom into account, one arrives 
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at a Schrodinger equation that can be solved by a separation ansatz, 



y, t) = ip^{x, t)^y{y, t) , (3) 
with (absorbing the separation constant into a state redefinition) 

■M.M = {^^ + ^-')M^,t), (4) 

We see from the first of these equations that ip* obeys a Schrodinger equation with standard 
kinetic term, but with the sign of the potential being reversed ('upside-down oscillator'). This will 
become important in the following discussion. Wave functions describing the Schwarzschild and 
the Reissner-Nordstrom black hole have been discussed in quantum geometrodynamics at various 
places, see for example pQ IHl UHl [HI [H] ■ This part of the total state is mimicked by the a;-system. 
Our purpose here is to construct Gaussian wave packets which are solutions of ^ and (O, and 
which describe the transition from the semiclassical regime to the final evaporation phase. For 
various relevant aspects of such wave packets in quantum mechanics, cf. [13 HH [E]. 
If ipoix, 0) denotes an initial state, the evolution for the x-part is found from 

dx' G{x, x'; t, 0)V'o(a;', 0) = ij{x, t) , (6) 

where G{x, x';t, 0) denotes the Green function. In the case of the inverted oscillator it reads (|18j. 
Sec. 6.2.1.8): 



G{x,x';t,0) = J -^^^Tf n 



impLUx- 



(x^ + x'^) cosh{LUxt) — 2xx' 



We shall now investigate the solutions of the Schrodinger equation for various initial states. 



(7) 



3.1 Squeezed Ground State 

Before we address the more interesting case of an initial coherent state in Sec. 3.2, we take in a 
first example as an initial state the ground state of the harmonic oscillator. 

For an ordinary oscillator, the system would stay in the ground state. Here, however, we find from 
dS]), together with the complex conjugation, cf. ([4]), the state 

rAx,t)^( ^, "^P""" ^'^ cxp (- (1 + i sinh 2uo,t)x^\ . (9) 

This is, in fact, a squeezed ground state. Comparing with the general form of Gaussian squeezed 
states, see for example [T7], one recognizes that the squeezing angle is = 7r/4, and the squeezing 
parameter is r = LVxt. The squeezing thus proceeds along the diagonal in phase space and increases 
linearly with time. Recalling the analogy to the black-hole case, this squeezing is expected to 
happen in the final evaporation phase when the kinetic term for the x-degree of freedom becomes 
significant. The squeezing vanishes in the formal limit mp — !■ oo. The evolution of the x-part of 
the wave packet is depicted in Figure 1. 

Assuming also for the y-part the ground state as initial condition, it is obvious that it will 
remain in this state if the standard harmonic oscillator propagator is used. The full state is thus 
the product of Q for the x-part with the standard ground state for the y-part. 
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The evaporation of a black hole would thus be mimicked by the squeezing of a Gaussian wave 
packet. It is known that Hawking radiation is in a squeezed state for semiclassical black holes 
[THmnj. Here, instead, we have a squeezing of the black-hole quantum state itself. 

This case of a squeezed ground state is unrealistic because the support of the wave function in 
the x-variable extends significantly into the negative regime. Therefore we shall turn to a coherent 
state as an appropriate initial condition. 

3.2 Initial Coherent State 

We shall now consider an initial coherent state for both the x- and the y-part. This choice 
captures more the idea of an initial semiclassical black hole in which the dynamical variable x may 
be assumed to be concentrated around an initial mass AI. Addressing first the x-part, we have 



where 




with Xq and po denoting the expectation values of x and p^, respectively. The state pop thus 
contains two free parameters xq and po: which together form the complex variable a. Application 
of ([7]) and complex conjugation then gives 

/. , hmpuJx smhujj - icoshw^jt |ap a*^ 1 \ 

exp la W — X — (- + icoshuJxt) \ . (12) 

\ V ft cosh2tJa;i 2 cosh2w^t^2 W ^ ^ 

Comparing p2p with ^ , we immediately recognize that this state experiences the same degree of 
squeezing. The state (jH) is, of course, a special case of ([T2l) from which it follows for a = 0. The 
absolute square of (fT2|) is depicted in Figure 2. 

Our model allows both signs for pq. In order to capture the idea of an initially decreasing 
mass we choose as initial condition a negative value for pq. This is shown in Figure 2. One can 
recognize that the localization of the packet centre moves toward smaller values of x before it 
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spreads. The spreading in this model should reflect the quantum gravitational behaviour in the 
realistic situation. In our simplified model it will be possible that also negative values of x can be 
reached, which is not possible if x represents the Schwarzschild radius; in a more realistic model, 
a potential wall at a; = should be introduced. 




The absolute square of (|12p can also be written in the form 

^^P ~^ Tt; — -{x - [xo coshuj^t J) , (13) 

\ ncoshZuJxt mpLUx J 

where the explicit form of f{xo,Po, t) is of less interest. It is easily seen from this result that the 
packet is peaked around the classical solution, but highly squeezed. We note that this fact makes 
it very sensitive to decohcrence, see Section 4. 

Taking for the y-system (the ordinary oscillator) an initial coherent state as in (|10p . one obtains 
the standard result for the time-dependent coherent state, as found in textbooks on quantum 
mechanics. 




The time evolution of the absolute square is shown in Figure 3. (More realistically, one should 
take into account a mild squeezing for this state, since it describes Hawking radiation [H].) The 
total state is then again a product of (|12|) and ([T4|) . 

3.3 Back reaction between a black hole and Hawking radiation 

In the previous section, back reaction was only implicitly implemented in the sense that the kinetic 
term was assumed to become relevant for small black holes. In order to be more realistic, a direct 
coupling between the black hole and its Hawking radiation, that is, a coupling between the x-and 
y-part, has to be included. This will mimic the process of back reaction more clearly. The simplest 
way is to include into ([2|) a linear coupling of the form fixy, where the parameter fi should depend 
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y 



Figure 3: Evolution of |V'5^{j/i under the ordinary oscillator propagator, with 
for simplicity. 



yo = POy 



on the black- hole parameters. Then equation ([2]), when considering only the a;- and y-part, will 
transform into 



iS^*(a;,y,t) 



92 



i^2mp dx^ 2my dy"^ 



(15) 



What would be a suitable initial state? As for the a:-part, one could still start with an initial 
coherent state. As for the y-part, mimicking Hawking radiation, we would suggest to take the 
following state as initial state, 



i'yoiV'to) oc exp 



2h 



-coth 



2TTL0yGM 

h lUJyto 



(16) 



where M is the original mass of the (Schwarzschild) black hole, which corresponds to the initial 
value xq of the x-part of the quantum state. What is the justification for this choice? As was 
shown in [7], p6|) corresponds to the state describing Hawking radiation. In [7] we have dealt with 
a dilaton model; for the Schwarzschild case one would expect the form (fTB|) of the state, cf. [TO] . 
It is appropriate to rewrite Equation (jlSp in the following way: 
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in which the following redefinitions are used: 

Qi cos 9 + Q2 sind 

X = 

y/mp cos 26 

Qi sin 9 + Q2 cos 9 



nl cos^ 29 



nl cos^ 29 



/niy cos 29 



uj~ cos 9 + iVy sni" 



2 • 2 

oj.,. sm I 



+ ujy cos^ ( 



/i sin 29 

jjL sin 29 
^mpuiy 



^iQi,Q2,t), 



Px = 



P, 



mp 
cos 26* 

cos 20 



Pi cos 9 + P2 sin 9 



Pi sin 9 + P2 cos 9 



cos2 29 
cos2 29' 



coth 



2'kujMM 



+ iwyio 



(17) 



(18) 
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Here, without any loss of generality, 



1 



cos2 20 \ 2 



^ {ojl+ujl)sin2e + ^ 



(19) 



must be zero in order that (fT7|) allows a separation of variables [21] • In this case we find 



/i = --^ympmy {ul + ujy) sin 20 



(20) 



with 6* g] — f,f[ in order to have the set of transformations physically consistent. It is 
important to point out that fi can assume any value in the interval 6 e] 
From (fT8|) and (|20p we obtain the following results. 



TT TT r 

4 ' 4 I 



1 



cos2 26* 
1 



iul i cos^ 6* - i sin^ 26 



^ cos2 29 



1 



sin^ 26* +ujI{ cos 



i sin^ 26* 



■ sin^ 29 



(21) 



We are now in a position to apply the same procedures as in Section 3. We must use, of course, 
the form of the initial states (fTU|) and (fTB|) under the redefinitions (fTH|) . In order to obtain the 
time evolution, the next step consists in calculating the following integral. 



J JdQ[dQ'^ G{QuQ'i;t)-G{Q2,Q'2;t) 



exp ( -^{Q'l sin^ ^ + Q'i cos^ 9 + Q[Q'^ sin 26*) 



1/4 



exp ( - (Qf cos^ 6* + sin^ 9 + Q'^Q'^ sin 26*) 



2u^ 

h 



{Q[ cos 6* + sin ( 



(22) 



where G{Qi,Q'i;t) is the inverted oscillator propagator and G{Q2,Q2it) the standard one. After 
some calculation we obtain 



mi,Q2,t)=[^) [-^J exp 



n^2 



2 W2 -^4 

2?i J^i 2h 



exp 



exp 



exp 



exp 



\2lox fi^iQicos9 fl2Q2siiit 
-ia ■ ' ' 



^1 



1^2(52 sin26' / l^iQi 



^3 



2.^1.^3 



ia* W — sinh flit cos | {ojy + ujx) 



ftsin 26'sinf]2i .^iQi , , /2tJx • , o ^ zil ^~ , ~ n2 
I — : ha y —— smh Sill cos y (w.y+ti^x) 



*2~ 



h \ h 

cos^ sinh r^it sin^ sin $12^ 



.^3 



jar 
2 



(23) 



8 



where we have defined the following functions of time, 

T\ = — ifii cosh ilit + Zjx cos^ Q sinh + u)y sin^ 9 sinh fJit, 

J-2 ~ — ^\ sinh r^ii — i^liujx cos^ cosh — if^icDj, sin^ cosh ^lit, 

JF3 = — iil2 cos Q.2t + i'l; sin^ sin + i^y cos^ 9 sin $12^, 

/— ~ \2 2 „ sinh r^it sin 1^2^ 

-(^jy+w^;) sm 26* — 

4y-i 

JF4 sin r22^ ~ zri2W2^ sin^ 9 cos $72^ — iVl2ijJy cos^ 9 cos 5^2^ 

„ ,~ ~ ^9 9 „ sinh rjit cos riot 

+ in2{ujy + LOxf sm^ 26* ^. 

47"! 



(24) 



Evaluating the functions J-"„ at t = enables us to recover the initial state ip{Qi, Q2, t = 0). Now 
we must perform the inverse coordinate transformation, 



J <3i = y/mpx cos 9 — y 
I (32 = ^frn^V cos 9 - ^ 

in order to give the explicit form of tp{x, y, t). Writing 



'myU sm f 
mpx sin ( 



(25) 



ipix, y, t) = F{t) cxp A{t)x'' + B{t)x + C{t)y'' + D{x, t)y 



(26) 



we find 



F{t) 



D{x,t) 



exp 



*2~ 



1/2 
-) X 

cos^ 6* sinh ilif sin^0sinri2i 



^1 



jar 
2 



exp 



a Uxiyiy + Wa;)^ cos^ ^sin 26* sinh sinri2t 



A(t) = - 



S(t) = 



^2 2 , 
— cos ( 



'4 . 2 

t;:^ sm 

^3 



{ujy + LOx) sin^ 20 



Of (wy +1^2,)^ sin^ 20 cos^ sin 1^2 i 



2'mpUJx 



ill cos^ 6* 5l2sin20 r22(i'y + ti3a;) sin^ 26* sinh flit 



4^1^3 



ill (wy + w^:)^ sin^ 20 cos^ 9 sinh flit sin Vl2t 



4^l2^3 

•^2 . 9 ^ -^4 9 ^ 17 1^2 

sin^ + T;r^ cos^ 9 + ^ 
^1 ^3 2^1^; 



■(tjj, + tiix) sin^ 20 



^\{ijJy + Wx)^ sin^ 20 sin^ sin 112^ 



X sin 20 



-y^mpuiy 



4^f^3 

——{ujy+ujx)cos29+ — + — 



Vi\ {ujy + cj^)^ sin^ 20 sin 112^ 



8^?^3 



ia* 2myUJx 



sin 20 



n2 



- (^y + ^x) COS 0sinhOi^- — + — 



ill (t^y + W:r)^ siu^ 20 siuh fli^ siu n2t 



(27) 
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The full state ((HI) describes an entangled state between this quantum black hole and its 
Hawking radiation toy model. One can then calculate from it the reduced density matrix for the 
black hole by integrating out the y-part. Tracing out the Hawking radiation thus gives a mixed 
state for the analogue of the black hole itself [7]. In fact, the presence of entanglement between the 
black hole state and the Hawking radiation is of particular interest as emphasized, for example, 
in [50]. It is important to recall that the black hole is an open system, which by itself (without 
taking into account the Hawking radiation as well as all other fields interacting with it) does not 
evolve unitarily. The information-loss problem for black holes can only refer to the closed system 
of black hole plus all other degrees of freedom which are entangled with it. By our very ansatz, 
in our model the full evolution of black hole plus Hawking radiation is unitary with respect to 
semiclassical time. 

The diagonal element of the reduced density matrix for the black hole is computed from 

Pxx^tT^yP = J \{x,y\x,y)\'^dy, (28) 

where we use \x,y) = ip{x,y,t). Inserting (^5]) into (^5]) . one finds, 

= \F\' cxp [x^A + A*) + x[B + B*)^ ^-^^exp (^^) ■ (29) 

Similarly, the off-diagonal clement of the reduced density matrix for the black hole becomes 

P..' = \F? exp [x'A + x'^A* + xB + x'B*^ ^-^^cxp , (30) 

where the sum D(x) + D*{x') implements an additional dependence on x and x' as we can verify 
from (|27p . Here we shall not discuss further the non-diagonal elements, which describe decoherence 
effects 0. 

Figure 4 depicts pxx for /i = 0, thus representing a situation where there is no explicit back- 
reaction term. The parameters are specified in order to have mp 3> my] this seems to be a 
reasonable choice, since the energy scale of Hawking radiation is much smaller than the Planck 
mass, except perhaps in the very last stage of the evaporation. We can observe that Figure 4 
is essentially identical with Figure 2 (obtained without taking into acount f{xo,po,t)). In fact, 
having /i = corresponds to take 9 = in ((27|) and ((29|) : consequently, we obtain the solution 
(ini), since lim = \il:^{x,t)\^ . 

9 — >0 

The most interesting situation, however, is to allow /i 7^ in order to have an idea of the 
effect of back reaction on pxx- Figure 5 shows what happens for ^ varying from to 100. In fact, 
significant modifications to pxx emerge mostly when ^ becomes larger than 1: we can observe a 
strong modification of pxx when the back reaction term comes into play. Due to the entanglement 
between the black hole state and the Hawking radiation, pxx is now sensitive to modifications of 
Wy or my. The particular choices of these parameters, which constitute the information carried 
by Hawking radiation, lead to very different forms of pxx- 

One recognizes from Figure 5 that there is for large /x, that is, for large back reaction, a 
supression of the strong squeezing. There remains a relatively narrow wave packet whose width, 
however, strongly oscillates. 

Instead of tracing out Hawking radiation from the full entangled state ([^S]). we can trace out the 
black-hole state and thereby arrive at the density matrix pyy. The computation of pyy processes 
in a similar fashion as for pxx, and the formal expression for it is comparable to ([59]) . The interest 
in this quantity arises from the fact that the effect of the back reaction between the black hole 
and the Hawking radiation (in the context of this simplified model) can be explored from what 
effectively leaves the hole. In particular, the issue of information being carried by this radiation 
is of special interest. The question we may ask is what should we expect to detect if the Hawking 
radiation were taken to be a signal affected by back reaction effects. In Figure 6 we represent pyy 
for various values of p. We recognize that the diagrammes for pxx and pyy look similar for large 
p. This could be due to the large entanglement between the black hole and its Hawking radiation, 
which leads to similar reduced density matrices. 
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Figure 4: Evolution of pxx, with mp = h = uJx = = 1 and po = —1; = for simplicity, and ^ = (no 
back-reaction coupling), ujy = ojx x 10^''^, rriy = nip X 10^^. In the contour plot the brighter areas correspond to 
higher values for pxx- 

4 Summary and conclusion 

The main purpose of this paper is to get an intuitive, though tentative, insight into how the 
final phase of blaek-hole evaporation may look like. In the geometrodynamical framework used 
here, the question is: suppose we have solved the quantum constraints and found the wave function 
describing the black hole and its radiation - which characteristics would identify it if some quantum 
state were detected? 

We have argued that in addressing this question we have to deal not with the full quantum 
constraints of the Universe, but with a (functional) Schrodinger equation which contains the exact 
Hamiltonian of the quantum black hole together with the semiclassical time t of the rest of the 
Universe. After all, there are observers who observe the quantum black hole from outside and who 
have a clock at their disposal. 

We have constructed and investigated a very simple model: two coupled harmonic oscillators 
which mimic black hole and Hawking radiation. We have found that our analogue of a black-hole 
state experiences a strong squeezing during evaporation, but that this squeezing may disappear 
for large back reaction; in this limit the reduced density matrices for black hole and Hawking 
radiation look similar to each other - their differences begin to disappear. 

The simplicity of our model allows at best to give some heuristic insight. It is most unlikely 
that the end phase of a black hole can be described by a system as simple as coupled oscillators. 
However, it is possible that some of the aspects discussed here - high squeezing of the black- 
hole state and large entanglement with Hawking radiation - occur also for the realistic quantum 
states. Moreover, features of an oscillator system are often discussed even in full approaches to 
quantum gravity, such as the recent proposal for a flux-area operator with an equidistant eigenvalue 
spectrum [5^. 

A big open question is the relation of such simplified models like ours to full field-theoretic 
approaches. This question can. of course, only be answered from the full theory, not from the angle 
of our model. But there exists an abundance of papers in the literature which address quantum 
aspects of black holes from this 'minisuperspace' point of view, with [51 [TUl EJ [T^ US] fisting only 
some of them. The hope there is always that one or the other result will survive in the full elusive 
description. Our paper is written in the same spirit as these papers. 

Many more questions remain to be answered. One concerns the calculation of the black-hole 
entropy and the recovery of the Bekenstein-Hawking formula in the limit of large black holes. Great 
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progress has been achieved here in loop quantum gravity [U [2] and string theory [3] , but some 
results have also been obtained from the Wheeler-DeWitt equation in quantum geometrodynamics, 
cf. [21] and the references therein. It would be of interest to give some exact results in the context 
of our simplified model. One must not forget that a black hole is a genuine open quantum system 
because it is susceptible to even small interactions with other fields (including its own Hawking 
radiation, as discussed here) [ID] . The process of decoherence would thus play a crucial role in the 
discussion [25j . This is especially important because it is known from quantum mechanics that 
squeezed states - and the final black-hole state in our model approaches such a state for not too 
large back reaction ~ are highly sensitive to decoherence. A possible connection with observation 
could be made when addressing primordial black holes - small relics from the early Universe, 
which in the appropriate mass range could be evaporating in the present phase of the Universe 
[26| . We hope to address some of these issues in future publications. 
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